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The paper reports on the quasi-static steady flow of a dry liquid foam around a fixed spherical 
bead, few times larger than the typical bubble size. The force exerted on the bead is recorded 
with a precision and a time resolution large enough to show the succession of elastic loading of the 
foam, separated by sudden force drops. The foam structure is observed by direct light transmission, 
synchronized with the force measurement, thus allowing to correlate the plastic events with the 
force variations. Scaling laws for the force signal as a function of the bubble size are detailed and 
interpreted with a simple elasto-plastic model. The spatial distribution of the plasticity is strongly 
localized in the first bubbles layers around the bead and the average size of bubble rearrangements 
increases with the corresponding force jump amplitude. 

PACS numbers: 83.80.Iz,83.60.La,62.20.Fe 



* Electronic address: Iisabelle.cantat@univ-rennesl.fr1 



2 



I. INTRODUCTION 

At small deformations, a dry liquid foam behaves like an incompressible soft elastic solid. The 
bubbles elongate and the amount of gas/fluid interface increases : the foam stores elastic energy 
and develops a restoring force. At larger strain, the deformed structure becomes instable, the bub- 
bles slide over each other and thus recover a more isotropic shape (for a review see Among 
the complex fluids, the foam is probably the material for which the elementary processes leading 
to elasto-plasticity are the best understood, and the numerical simulations the more accurate. The 
involved scale is the millimetric bubble scale, at which the structure and/or the velocity field can 
be monitored with various resolutions and various techniques, among which X-ray tomography]^, 
diffusing wave spectroscopy (DWS) optical tomography 6, 7,^J, magnetic resonance imaging 
(MRI) la [13, Liquid foams are thus very promising model systems to understand the cou- 
pling between the local organization and the macroscopic rheological response. Generalization is 
expected to other complex fluids, as colloidal suspensions, slurries, pastes or clays, in which the 
structure scale is nanometric. 

Most of the experimental results on 3D foams concern homogeneous simple shear geometries. 
The intermittency, the stress and strain inhomogeneities and their tensorial nature are crucial 
parameters that cannot be tackle with classical rheometers. In contrast, the Stokes experiment, ie 
the motion of a rigid obstacle in a liquid system, induces an inhomogeneous stress in the material 
and thus provides very complementary information. Some previous works report on the relation 
between the velocity of a foam flow and the force it exerts on an obstacle in 2D or in 3D 

15.J 16]. Nevertheless, the force fluctuations, correlated to the direct observation of the rapid 
bubbles rearrangements responsible for the plastic transformations are presented here for the flrst 
time. 

Our experimental set-up, detailed in section^] imposes a steady foam flow around a bead. The 
force exerted on the bead is recorded with a precision scale under which the bead is rigidly fixed, 
as already used in . A camera is synchronized with the scale and records the images of the foam 
structure. In section UTTI the foam shear modulus and its plastic threshold are deduced from the 
force signal and compared to predictions obtained with a continuous elasto-plastic medium. The 
force fluctuations are detailed in section Hvl and the correlation with the magnitude of the plastic 
events, obtained from the image processing, is demonstrated. Finally, the spatial distribution of 
the foam structure rearrangements is discussed in section fVl 

II. EXPERIMENTAL MEASUREMENT 

A monodisperse dry foam, with millimetric bubbles, is produced in situ in a transparent 10 x 
10 X 20 cm Plexiglas box by nitrogen blowing in a solution of SDS (sodium dodecyl sulfate) 3g/L 
and dodecanol O.Olg/L. The bubble size can be varied by changing the gas flux and the injector 
needle diameter. The surface tension of the solution is 7 = 37mN/m, its viscosity rj — 10^^ kg m^^ 
s~^ and its density p = 10'^ kg/m'^. Once the cell is filled, the box is closed to avoid evaporation 
and films breakage, and the remaining solution is removed from the cell bottom. As discussed in 
|16|. aging processes (drainage and coarsening) strongly perturb the measurements so we let the 
foam draining during half an hour, until the liquid fraction becomes constant and of the order 
of 1 or 2 %. The coarsening is slowed down by use of CeFi4 traces added to the nitrogen. The 
experiments presented in this paper have been done between half an hour and two hours after the 
foam formation in section Hill and between half an hour and one hour in section Hvl 

All the bubbles have the same volume, with a dispersion smaller than 10%. Its value, which is 
our main control parameter, can be obtained by taking a picture of a spherical bubble rising in 
the liquid phase, during the foaming process. Nevertheless, it is more convenient to measure the 
average distance d between two horizontal films within the bubble column at the box corner. The 
first method has thus been used only once to determine the relation between d and the bubble 
volume, or, as used classically in the literature, the radius ry of the sphere of same volume. We 
found 



rv = 0.54 d± 5%, 
which will allow quantitative comparisons with previous works. 



(1) 
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A glass bead of radius R = 0.5cm is maintained at a fixed position, in the middle of the horizontal 
box cross-section, by a rigid glass rod. The sphere and the rod are entirely wetted by the solution. 
The box can translate vertically at a very low velocity V and thus imposes the foam motion around 
the immobile bead (see Fig^. The rod is fixed below a precision balance (Sartorius BP 211D) 
and the force exerted by the foam on the bead is recorded every 0.25s with a precision of 1 /iN. 
The reference force is taken without foam in the box. 

When the box begins to move, the force firstly increases during a transient, which is disregarded 
in this study (see Fig|2Jl. The end of this transient is estimated roughly as the time at which 
the force reaches its mean value and the data collected before this time are suppressed from 
the statistics. Then we observe a succession of sharp force decreases (plastic events) separated 
by smooth force increases (elastic loading) during which the force varies linearly with the box 
displacement x, even for the longest force increases observed (see FigEJ. The beginning and 
the end of the elastic loading indexed by i are denoted respectively by tmin,i and tmax,i and are 
extracted from the raw force data with the following method. After the jump i — 1, a new loading 
i begins at the first time verifying F{t + dt) > F{t) + e, with dt=0.25s. F{tmin,i) is denoted 
by Fmin,i- The parameter e is fixed slightly above the noise level at e = 5/iN. The end of the 
loading i (or equivalently the beginning of the jump i) occurs at the first time tmaxA at which 
F{t + dt) < F{t) — e, and Fmax,i = F{tmax,i)- If the foam moves upwards, we replace F by -F 
in the data processing. From this list of force maximum and minimum, we compute the jump 
amplitude 5Fi = F„iax,i-i — Fmin,i, the rate of force increase during elastic loading dF/dt)i = 

{Fmax^i - Fjmn,i)/{t7nax,i " tjmn,i), the Waiting time between two jumps t^^i = t^ax,i - tmax,i~l 

and the jump duration td^i = tmin,i — tmax,i-i- The time resolution is much smaller than the 
typical force jump duration (of the order of 1.5 s as demonstrated in section llV|) and is not a 
limiting factor. The parameter e has been chosen in order to disregard the small amplitude force 
jumps for which the correlation with the plastic events observed on the images disappears. We 
checked carefully that it does not influence the amplitude distribution of force jumps larger than 
e. The end of the transient is not precisely determined, but the uncertainty thus induced on the 
mean force value is smaller than the error bar and it does not modify the other quantities. 

As demonstrated on Figs. Eland[7| the force signal is symmetric when the displacement direction 
is modified and the presence of the rod seems to be without major infiuence. The Archimedes 
force, that could also break the up/down symmetry, is of the order of 4 / 3 n g p(j> 50/iN, which 
is much lower than the typical force measured, of the order of few mN. The flow occurs in a 
quasi-static regime, and, as checked experimentally over the range (1-50/im/s), the velocity only 
rescales the time unit. This is consistent with the small viscous force value estimated as follow. 
The viscous forces /„ between the foam and the bead are localized along the Plateau borders and 
are of the order of iOjCa^^^ ^ 20/iN/m per unit length of Plateau border touching the sphere, 
with Ca = r]V/'f ^ 10~^ the capillary number [T7|. The resulting force is then approximately 
SwR^fy/d ^ 3/iN, which is negligible. Instead of dF/dt, the pertinent force variation is thus 
dF/dx, with x the box displacement, simply obtained using the relation dF/dt = VdF/dx. 

The foam is lit by a cold and uniform light of the same size than the box (Schott-Fostec DDL) and 
the foam images are recorded from the opposite side with a camera focused on the bead, at a rate of 
4 images per second. The Plateau borders diffuse the light and appear in black on white, projected 
on a 2D image of the whole 3D foam (see FiglSJ. The slow motion of the foam at velocity 5/im/s 
induces a translation of 1.25/im between two successive images, which is much smaller than the pixel 
resolution {QQ^rn) and than the Plateau border typical radius (a = 200 — 300 /xm). Therefore, if the 
foam only translates in front of the camera, the difference between two successive images shows very 
low grey levels, corresponding to the noise. In contrast, bubbles rearrangements by Tl processes, 
ie neighbor exchange between bubbles, induce rapid fllms motions and lead to the typical patterns 
shown on Fig^Ja). For bubbles larger than 2.5 mm, these plastic transformations can be localized 
automatically by the following computer image processing, made with the commercial software 
Visilog. The difference between two successive images is binarised, keeping only the largest and 
smallest grey levels, smoothed on the bubble scale and binarised again. A succession of dilatation 
and erosion allow to keep only domains with a size comparable to one bubble size or larger. We 
finally obtain domains without hole called " spots" in the following. They are approximately of the 
size of the few bubbles involved in the corresponding plastic event (see Fig0Jb)). 

The typical duration of a plastic event is of the order the second and it is detected in average 
on 1 or 2 images. The criterion to identify two successive spots as the same event is that their 
centers of mass are closer than the geometric average of their characteristic radius, computed from 
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their areas. Because we have only a 2D information, it may happen that two distinct events lead 
to spots at the same place on an image or on two successive images. Nevertheless, as discussed 
in section IIVI only 2 or 3 plastic events happen during a force jump so this spots superposition 
occurs rarely and the case has been neglected. The presence of a plastic event at a given time is 
well detected by our method, except for very small displacements that usually correspond to the 
beginning or the end of a larger event. In contrast, it is difficult to define and to quantify precisely 
the event size, which depends obviously on the grey level threshold and of the pixel resolution. The 
slightly arbitrary volume assigned to a plastic event in the data presented below is the area of the 
corresponding spot to the power 3/2 (or, in case of an event with successive spots, the area of the 
largest one). It can be used to compare plastic events between each other, but not to determine 
the number of elementary Tl events involved (each Tl involving only four bubbles). For further 
experiments, the use of two synchronized crossed cameras will improve the determination of the 
size and position of these rapid rearrangements. Finally, no information is obtained from the foam 
column in front or behind the bead. 



III. SHEAR MODULUS AND PLASTIC THRESHOLD 

Among other technics summarized in tabled the foam shear modulus is usually measured in a 
simple shear geometry, by small amplitude oscillations around an equilibrium position, as the ratio 
between the imposed displacement and the macroscopic stress. In contrast, the Stokes geometry 
induces non-uniform stress and strain in the foam, in a steady state flow. Nevertheless, assuming 
a purely elastic response of the medium during the smooth force increase, an expression for the 
shear modulus can be extracted from the comparison between the measured force and a continuum 
medium analytical predictions. 

The distribution of the slopes of the force-displacement data during the elastic loading is rel- 
atively narrow (see Fig. |21and[SJ). The mean slope was computed for each flow using a weight 
proportional to the duration of each slope (or, equivalently, all individual force increases are added 
without weight and the sum is divided by the total time spent in the elastic state). The exper- 
imental data are in good agreement with the expected linear dependence with l/d, for a bubble 
size varied by more than a factor 2, and is well adjusted by the relation given on FigEl 

dF _ 2.0 10-3 ± 15% 

'dx ~ d ^ ^ 

with the force slope expressed in N/m and d in m. The error bar obtained from the best Gaussian 
fit displayed on Fig. [S] and the incertitude on the fitting parameter on Fig. |B| both lead to a 
precision of the order of ±15% 

The force needed to displace a rigid sphere embedded in an infinite incompressible continuum 
medium is simply, assuming a free sliding on the bead, 

Felas ^^TtGRS^ , (3) 

with R the sphere radius, Sx the displacement amplitude and G the shear modulus of the sur- 
rounding medium (see appendix A). The foam shear modulus is thus deduced from eqs[21and|3 
Once scaled by j/d or j/ry (see eq^J, we obtain : 

G = 0.86 7/d ± 15% = 0.46 7/rv ± 20% . (4) 

This value is compared to previously obtained results in the table Slip at the solid boundaries 
is often mentioned to explain the dispersion of the experimental results obtained with classical 
rheometers. A better control is obtained with our set-up that assumes a perfectly slipping surface 
on the bead, which is easier to ensure than a no-slip condition. Nevertheless, the discrepancy 
between the various results may arise from more fundamental reasons related to geometrical effects. 
Our value is indeed very close to the value obtained by the indentation test Ts'l , which is also based 
on a local mechanical action on the foam. Numerical simulations on the Kelvin cell done by Kraynik 
et al demonstrate that the foam shear modulus varies between 0.357/ry and 0.607/rv, depending 
on the shear orientation p^ . Structure reorganizations around the bead, due to the many Tl 
induced by the flow, may lead to an important local foam anisotropy and to a decrease of the shear 
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modulus in the direction of the applied force. 

The yield stress is more difHcult to define than the shear modulus. We deduced it from the mean 
force value, during the steady state. This force Fth is found to varies linearly with 1/d (see Fig. 
[Tj) and is well fitted by the relation, scaled by Rl^j/d, 

Fth = ^■65R'^"f/d±15% . (5) 

In the frame of the purely elastic model presented above, the stress tensor has a simple analytical 
expression, detailed in Appendix A. The shear stress, which can be more easily compared to 
previous values obtained in classical rheometers, reaches its maximal value amax in front of the 
bead or at the rear : 

amax — 1-5 — - — 0.75 — . (6) 

From eqsI21E]and|ll we deduce the mean value of the effective bead displacement in the steady state 
(with respect to its equilibrium position, deduced from the purely elastic model) Sx,th = 0.52R. 
Using eqlHl we finally get the yield shear stress, scaled by j/d, 

= 0.65 7/d ± 15% = 0.35 7/rv ± 20% . (7) 

Previous experimental results are very scattered for dry foams, with a measured yield stress in the 
range (0.04 - 0.2 "f/rv) (see 113,01 and references therein). Gardiner et al also mention the value 
^th = 0.56 7/ry for a very dry foam '2^1 . Numerical simulations by Kraynik et al on polydisperse 
foams lead to Utn = 0.187/ry with = 0. As discussed in the results scattering is probably 
mainly due to a rapid variation of uth with the liquid fraction when 4> is close to and to the 
various definitions used, which have not been demonstrated to be equivalent. The surprisingly 
large value we obtain may be due to the fact that the maximal stress is only reached in a very 
small foam volume, which decreases the probability to found a weak structure in the concerned 
volume. Numerical simulations of a Bingham fluid or an Herschel-Bulkley fluid flowing around a 
sphere also leads to a relation between ath and f|2^|23,|3|- Anyway, as a non-slip condition is 
assumed on the sphere, a direct comparison is not possible. 

The maximal elongation of the first bubble layer before detaching from the sphere is another 
characteristic of the plastic threshold. Unfortunately the image resolution is not sufficient to obtain 
a direct measure. Within the continuum elastic medium approach, the maximal elongationnal 
stress Grr = '2-G5x/R is determined in Appendix A, using the spherical coordinates with a polar 
axis oriented in the direction of the bead motion. The maximal elongationnal strain e^r at the rear 
of the bead, is thus : 

e„ = Urr/i^G) = 0.52 ± 10% , (8) 

The limit of stability of a 2D foam between two plat es when the gap width is increased has been 
recently determined analytically by S. Cox et al [23. A n-sided bubble detaches from one of the 
two plates when the gap g between the plates reaches g{n) = 1.1 ry\/6/(6 — n), leading to the 
maximal elongation for n = 5 of the order of (^(5) — 2ry)/(2rv') ^ 0.35, of the same order of 
magnitude as our experimental result. 



IV. FORCE JUMP AND Tl EVENTS. 



The elastic loading is interrupted by sudden force decreases. For young foams, the force jumps 
distribution is well fitted by an exponential law (see Fig. (S)). This result questions the universality 
of the power-law behavior of the stress drop found in 2D Couette geometry. A power law with 
exponent -0.8, followe d by an exponential cut off, was found experimentally by Dennin et al with 
a sheared bubble raft |26l l27l | , in good agreement with numerical results obtained by Durian et 
al using the bubble model (exponent -0.7) 123. 129L l30l |. Okuzono et al, using a vertex model, also 
found a power law, but with a exponent -1.5 |3l| . In our experiment, the force jump distribution is 
qualitatively modified when the coarsening process becomes non-negligible and a power law with 
an exponent close to -1.5 is then observed (see Fig. ISJ. The detailed study of the transition between 
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both behaviors will be the aim of a future work. This paper only focuses on the non-coarsening 
foams and the oldest foams, that however have the same shear modulus and plastic threshold than 
the youngest, are disregarded in the following. 

If the mean force jump < SF > is simply computed from the list of recorded force jumps, its 
value depends on the arbitrary parameter e, which is the cut off for the small jumps. We thus 
chose to extract < SF > from the exponential fit, extrapolated to zero. We checked that both 
methods lead to results variations smaller than the error bar. The averaged jump amplitude as a 
function of the bubble size is well fitted by the linear relation (see Fig. llUfl : 

< 6F >== l.Ojd ± 10%. (9) 

This scaling law is consistent with the simple approach which identifies the plastic event responsible 
for the force jump with a single bubble detaching from the bead. Neglecting the contribution of 
the pressure (which may in fact play a significant role, at least in 2D 0|), the force exerted on 
the bead roughly scales as — C~ , with £+ and C~ the total length of Plateau borders touching 
the sphere respectively on the downstream side and on the upstream side. As detailed on Fig. 
these lengths scale as ~ RVu^, with the numbers of bubbles on each side. If one bubble 
detaches from the downstream side, the variation of scales as SC^ ~ Vn+ ^ d which leads 
to 5F ^ 'fd. The foam behavior can thus be understood, at the first order, with very localized 
events, which is confirmed by the spatial Tl distribution discussed in the section Ivl 

The mean force jump duration is < tmax,i ~ tmin,i+i >= l-5s ± 0.5s. A closer analysis of the 
force signal indicates that the relevant criterion for the jump beginning is the slope modification 
and not the maximum of the force. This is quite difficult to detect automatically and thus an 
arbitrary duration has been assigned to the force jump i, which begins at tmax,i — 2(s) and finishes 
at tmin,i+i + 0.5(s). With this convention a satisfying pairing has been obtained between the 
images and the force measurements, for all the experiments made with bubbles larger than 3mm 
in diameter (for smaller bubbles no information is obtained from the images). In average, 85% 
of the spots occur during the 15% of the time corresponding to force jumps. The 15% remaining 
spots probably correspond to small force jumps, smaller than 5^N, of the order of the noise and 
not taken into account. We obtain in average 2.5 plastic events per force jump. For each force 
jump, the volumes of these events, determined as explain in section^ are summed to get the total 
volume of foam reorganized during the jump. The Fig ll II shows the correlation between the force 
amplitude and the volume of reorganization. 

Most of the plastic events only appear on one image, which makes the estimation of the event 
duration difficult. Nevertheless, the number of successive images on which a given plastic event 
appears is exponentially distributed, and the average event duration deduced from this distribution 
is T = 0.2s. This time is much smaller than the time scale obtained from the force jump. The 
Plateau borders appear on the image with a typical size of 5 pixels and it moves from its own 
thickness between 2 images only if the velocity is of the order of 2 mm/s. This large velocity is 
probably only reached during a small fraction of the whole Tl transformation, thus explaining the 
difference between the two time scales. 



V. PLASTIC EVENTS SPATIAL DISTRIBUTION 

The correlation between the localization of the plastic events and the force jumps amplitudes 
would be very interesting, and will be feasible with a better jump statistics and an improved 
determination of the spot position in the 3D space. In this paper, we simply discuss the spatial 
distribution of the plasticity, averaged over the whole measurement time. This does not require 
to identify each individual event and thus allows to suppress the thresholding from the image 
processing and to take advantage of the whole information contained in the grey levels. After 
subtracting two successive images, we take the absolute value and we put to zero all the pixels with 
a grey level smaller than the noise level. The other grey level values, containing the information 
on the displacements, possibly smaller than one pixel, remain unchanged. Then all the images are 
summed. After a grey level binning process keeping only 5 different values, only used for a better 
legibility, we obtain typically the image shown on Figll2l 

The isodensity lines are organized around the two high stress regions, in front of the bead and 
behind. The experimental data correspond to a projection in the (x,y) plane, but, as the plasticity 
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decreases quickly with the distance from the bead center, it has quaUtatively the same angular 
behavior as the distribution in the median plane defined as z — (the frame orientation is shown 
on Fig. I12|l . It can thus be directly compared to the maximal shear stress angular distribution 
in the z — plane for a purely elastic medium, as plotted on Fig 1121 Both distributions show 
similarities, but the plasticity lobe upstream {x > 0) extends further laterally than the other 
one (independently of the presence of the rod), whereas the elastic shear stress distribution is 
symmetrical (see Fig ll2|) . 

The radial distribution has been investigated along the line y = which has the largest sig- 
nal/noise ratio. Local fluctuations have been reduced by averaging the grey level along a 25 
pixels segment oriented in the Uj^ direction. The grey level decreases phenomenologically as 
fci exp(— as depicted on Fig ll3l As the projection only induces polynomial transforma- 
tions of the signal, the exponential decreases also applies in the 3D space, along the line (y=z=0). 
In a purely elastic medium, the stress amplitude decreases as along the same line (see ea lT^ . 

The typical distance over which the plasticity disappears is of the order of 2d ± 0.3 and seems 
to be slightly larger in presence of the rod. The plastic events are thus almost all in the first 
and second layers of bubbles and the foam structure is not strongly modified further. Leonforte 
et al performed numerical simulations in a 3D Lennard-Jones fiuid in a similar geometry, in the 
linear regime 's^] . They demonstrate that the average stress field is well described by the classical 
elasticity in this case, but that the stress fluctuations decrease exponentially when the distance to 
the bead increases. This may be related to the exponential decrease we obtained experimentally, 
but the characteristic distances obtained numerically are much larger, of the order of 20 interatomic 
distances. Exponential decrease of the plastic events density has already been observed in 2D foams 
in Couette geometry, with a typical distance of the order of few bubble layers 

VI. CONCLUSION 

This article reports on a Stokes flow experiment performed in a 3D dry foam, in a quasi-static 
regime. We provide an original measure of the foam shear modulus and of the plastic threshold, 
and show that a simple model of elasto-plastic material allows to predict the scaling laws obtained 
for the mean force exerted on the bead, as well as for the average value of the force fluctuations. 
The structure reorganizations are recorded and their characteristic volumes are correlated to the 
amplitude of the force jumps they induce. Finally we show that the plastic events density decreases 
exponentially with the distance to the bead, as it is typically observed in case of shear banding in 
foam Is^l or in granular matter [s^. 

In conclusion, this experiment is an interesting alternative to the more classical simple shear 
geometry to probe the mechanical response of the foam. As the force is localized, only a small 
amount of bubbles are involved and the stress and strain fluctuations are more pronounced. Ad- 
ditionally, there is no need to suppress the sliding along the wall : the bead surface is smooth 
and thus induces no external random noise. The preliminary results concerning the influence the 
coarsening on the distribution of the force fluctuations and on the plasticity localization will be 
studied in greater details in a future work. 
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Appendix A 

For an elastic medium of Lame parameters A and G, the constitutive relation between the stress 
tensor aij and the deformation tensor e^j = l/2(diXj -\- djXi), with di the first derivative with 
respect to the i*'* spatial coordinate, is 



CTy = ATr(e) + 2Ge,j (10) 
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The equilibrium equation, without external force is given by 

= Gj^^jjXi + (A + G)d,{divX) 



(11) 



In the incompressible limit, X oo, it becomes 

= G ^ djjX, ~d,P , = divX 



(12) 



with P = -l/3rr(CT) = - (A + 2G/3)Tr(e) 

These equations are analogous to the hydrodynamical equations in the viscous regime. The 
displacement field around a sphere of radius R displaced from its equilibrium position on the 
distance 6x in an infinite 3D incompressible elastic medium is given, in spherical coordinates 
(r,0,,/)),by m 



Xr = Sr cos 9 — 



R 

Xe = -Sx sm0— 
2r 



(13) 



The direction u^; of the bead motion is the polar axis of the spherical coordinates and the bead 
center is the frame origin. A vanishing tangential force is assumed on the sphere. The resulting 
stress tensor is 



a = -PI + G5, 



RcosO 



-2 
1 
1 



The maximal shear stress value, relatively to the orientation of the normal vector, is 

Rcos9 



ash = 1.5 GSx- 



(14) 



(15) 



and the resulting force on the sphere is 



F = AttGRSx 



(16) 



[1] A. M. Kraynik, Ann. Rev. Fluid Mech. 20, 325 (1988). 

[2] D. Weaire and S. Hutzler, The physics of foams (University press, Oxford, 2000). 

[3] R. Hohler and S. Cohen-Addad, J. Phys.: Condens. Matter 17, R1041 (2005). 

[4] J. Lambert et al, Colloids and Surfaces A 263, 295 (2005). 

[5] D. J. Durian, D. A. Weitz, and D. J. Pine, Science 252, 686 (1991). 

[6] C. Monnereau and M. Vignes-Adler, Phys. Rev. Lett. 80, 5228 (1998). 

[7] F. Rouyer, S. Cohen-Addad, M. Vignes-Adler, and R. Hohler, Phys. Rev. E 67, 021405 (2003). 
[8] A. van Doornum and S. Hilgenfeldt, preprint (2006). 

[9] B. Prause, J. Glazier, S. Gravina, and C. Montemagno, Condens. Matter 7, L511 (1995). 
[10] P. Coussot et al, Phys. Rev. Lett. 88, 218301 (2002). 
[11] S. Rodts, J. C. Baudez, and P. Coussot, Europhys. Lett. 69, 636 (2005). 
[12] B. DoUet et al, Phys. Rev. E 71, 031403 (2005). 

[13] B. DoUet, M. Aubouy and F. Graner, Phys. Rev. Lett. 95, 168303 (2005). 

[14] S. J. Cox, M. D. Alonso, S. Hutzler, and D. Weaire, The Stokes experiment in a foam (MIT-Verlag, 

Bremen, 2000), pp. 282-289. 
[15] J. R. de Bruyn, Rheol. Acta 44, 150 (2004). 

[16] L Cantat and O. Pitois, J. Phys.: Condens. Matter 17, 3455 (2005). 
[17] L Cantat, N. Kern, and R. Delannay, Europhys. Lett. 65, 726 (2004). 

[18] M. F. Coughlin, E. P. Ingenito, and D. Stamenovic, J. Colloid Interface Sci. 181, 661 (1996). 
[19] A. M. Kraynik and D. A. Reinelt, J. Colloid Interface Sci. 181, 511 (1996). 

[20] B. S. Gardiner, B. Z. DIugogorski, G. J. Jameson, and R. P. Chhabra, Journal of Rheology 42, 1437 
(1998). 



9 



F. Rouyer, S. Cohen-Addad, and R. Hohler, Colloid. Surf. A 263, 111 (2005). 

A. N. Beris, J. A. Tsamopoulos, R. C. Armstrong, and R. A. Brown, J. Fluid Mech. 158, 219 (1985). 
M. Beaulne and E. Mitsoulis, J. Non-Newtonian Fluid Mech. 72, 55 (1997). 

J. Blackery and E. Mitsoulis, J. Non-Newtonian Fluid Mech. 70, 59 (1997). 

S. Cox, D. Weaire, and M. F. Vaz, Eur. Phys. J. E 7, 311 (2002). 

J. Lauridsen, M. Twardos, and M. Dennin, Phys. Rev. Lett. 89, 098303 (2002). 

E. Pratt and M. Dennin, Phys. Rev. E 67, 051402 (2003). 
D. J. Durian, Phys. Rev. Lett. 75, 4780 (1995). 

D. J. Durian, Phys. Rev. E 55, 1739 (1997). 

S. Tewari et al, Phys. Rev. E 60, 4385 (1999). 

T. Okuzono and K. Kawasaki, Phys. Rev. E 51, 1246 (1995). 

F. Leonforte, R. Boissiere, A. Tanguy, J. P. Wittmer, and J.-L. Barrat, |co nd-mat /05056io] (2005). 

G. Debregeas, H. Tabuteau, and J.-M. diMeglio, Phys. Rev. Lett. 87, 178305 (2001). 
O. Pouliquen and R. Gutfraind, Phys. Rev. E 53, 552 (1996). 

E. Guyon, J. -P. Hulin, , and L. Petit, Hydrodynamique physique (EDP science, Paris, 2001). 

B. Derjaguin, KoUoid Z. 64, 1 (1933). 

B. Derjaguin and E. Obuchov, KoUoid Z. 68, 16 (1934). 

D. Stamenovic and T. A. Wilson, J. Appl. Mech. 51, 229 (1984). 

H. M. Princen and A. Kiss, J. Colloid Interface Sci. 112, 427 (1986). 
B. Budiansky and E. Kimmel, J. Appl. Mech. 58, 289 (1991). 

D. Stamenovic, J. Colloid Interface Sci. 145, 255 (1991). 

Q. Sun, J. P. Butler, J. P. Suki, and D. Stamenovic, J. Colloid Interface Sci. 163, 269 (1994). 

T. G. Mason, J. Bibette, and D. A. Weitz, Phys. Rev. Lett. 75, 2051 (1995). 

D. A. Reinelt and A. M. Kraynik, J. Fluid Mech. 311, 327 (1996). 

A. Saint- Jalmes and D. J. Durian, Journal of Rheology 43, 1411 (1999). 



10 



TABLE I: Previous experimental and theoretical results for foams or emulsions shear modulus in the dry 
limit. 
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Derj aguin [37] 


E, foam 


rheometer osc. 
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Stanienovic[38] 
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rh. stat. 
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Sun [42] 


E, foam 


Wave prop. 
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rh. osc. 
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Coughlin[18] 


E, foam 


Indentation 


0.4 


KravnikflQ. 441 


N, Kelvin C. 


Surf. Evol. 


0.35-0.60 


Saint- Jalmes [451 


E, foam 


rh. osc. 


0.51 


C. and P. 


E, foam 


Stokes 


0.46 
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Figure captions 

Fig. 1 Experimental set-up. 

Fig. 2 Typical force graph obtained with the foam moving at 5 /im/s. The box is moving 
downwards, so the force exerted by the foam on the bead is oriented downwards too and the 
balance signal is positive. The transient finishes at t = 1000s. The open circles on the magnified 
graph represent the maxima and minima recorded and the dots are the measurement points. 

Fig. 3 Raw image obtained from the camera. The small rectangle is shown on Fig^ at two 
stages of the image processing. 

Fig. 4 A Tl event as it appears on the difference between two successive images (left), and 
after image processing (right). The value retained to measure the plastic event volume is the area 
of the spot to the power 3/2. 

Fig. 5 Probability distribution of the force variation with the bead displacement in a 3mm 
bubbles foam, during the elastic loading. It is obtained from the slopes of all the smooth force 
increases of a typical flow, weighted by the force increase duration. The events number is 162. 
The full line is the best Gaussian fit. 

Fig. 6 Mean value of the force slope as a function of the bubble size d. The simple dots and the 
dots with a circle correspond respectively to an upwards and to a downwards motion of the bead. 
No significant difference is observed between both directions of motion, so the asymmetry due to 
the rod is less than our error bar (of the order of 15%). The full line is the best fit of the type y=ax. 

Fig. 7 Force exerted by the foam on the bead, averaged over the steady state, (o) : downward 
motion of the box, (x) : upward motion. The solid line is a linear fit. 

Fig. 8 Force jumps distribution for a 45 min old foam with 2.8mm bubbles (no coarsening, 
68 events) and for a 2h old foam without CqFi4 (coarsening, 251 events) in semi-log and log-log 
(in insert) representation. The event density Pjump(^) related to the number of force jumps 
^*jump which have amplitude comprised between A — AA/2 and A + AA/2, observed during 
a 1 micron bead displacement, by p^^^-^p{A)AA = "jump- Only the force jumps larger than 
10 fiN are recorded. We used a non-constant binning to ensure that every point represents 
approximatively the same events number. The horizontal error bar is thus the distance between 
two points in the x-direction. A power law fit seems to be more appropriate for the coarsening 
foam (Pjump(^) ^ A^^-^^°-^/im^^mN^^ with A in mN), whereas an exponential decay well 
agrees for the non-coarsening foam (pjunip(^) ~ io^/0-2i±o o5 -^^jth A in mN). The error bars on 
the exponents are deduced from the extremal slopes compatible with the horizontal error bars on 
the data points. The first point is disregarded for the coarsening foam. 

Fig. 9 Schematical representation of the force variation induced by a plastic event, are 
the numbers of bubbles on each sphere side (upstream and downstream) and are the typical 
diameters of their contacts with the sphere. The relations n^{d^y ^ and rv^d^ ^ leads 
to the relation ^ rVu^ used in the text. 

Fig. 10 Mean jump force computed from an exponential fit of the force jump distribution, 
extrapolated to 0, in order to take into account the very small jumps, despite the noise. The data 
are in good agreement with a linear fit (full line). 

Fig. 11 Average volume of foam reorganized during a force jump as a function of the jump 
amplitude, for a 3mm bubble foam. The averages are performed after a force jumps binning and 
the number of jumps corresponding to each point of the graph is proportional to the circle radius 
(110 force jumps were recorded in total). For each force jump amplitude, the standard deviation 
of the distribution of reorganized volumes is depicted by the error bars. The insert shows the 
same data in a log-log plot. 
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Fig. 12 Spatial repartition, projected in 2D, of the foam quick displacements, averaged during 
a 1cm displacement of the box, upwards (— u^, direction). The averaging process is discussed in 
the text. Light and dark grey levels represent respectively the high and low density of plasticity. 
The full lines are the curves of constant maximal shear stress value in the median plane of a 3D 
elastic medium, obtained from eq. 1151 

Fig. 13 Radial variation of the plastic events density. Series C and D : The grey levels shown on 
Fig ll2l are averaged over 25 pixels on the lateral direction, along two vertical lines passing through 
the bead center, in front of and behind the bead (before the grey level binning). The resulting 
grey levels are plotted as a function of (x/d)'^, which underlines the good agreement with the 
phenomenological adjustment A:i exp(— a;^/fc|), with fci an arbitrary prefactor and k2 the typical 
distance over which the plastic events density decreases. The origin for x is the center of the sphere. 
Series A and B : same kind of data for a foam moving in the other direction, shifted vertically for 
a better legibility. 
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FIG. 3: Cantat, Physics of Fluids. 
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FIG. 5: Cantat, Physics of Fluids. 




FIG. 6: Cantat, Physics of Fluids. 
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FIG. 7: Cantat, Physics of Fluids. 



20 




Force jump (mN) 



FIG. 8: Cantat, Physics of Fluids. 
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FIG. 10: Cantat, Physics of Fluids. 
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FIG. 11: Cantat, Physics of Fluids. 
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